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Abstract. We consider a branching process with Poissonian immigration where individuals 
have inheritable types. At rate 9, new individuals singly enter the total population and start a 
new population which evolves like a supercritical, homogeneous, binary Crump-Mode- Jagers 
process: individuals have i.i.d. lifetimes durations (non necessarily exponential) during which 

t— I they give birth independently at constant rate 6. First, using spine decomposition, we relax 

previously known assumptions required for a.s. convergence of total population size. Then, 

we consider three models of structured populations: either all immigrants have a different 

fi type, or types are drawn in a discrete spectrum or in a continuous spectrum. In each model, 

the vector (Pi, P2, . . . ) of relative abundances of surviving families converges a.s. In the first 

f~^\ model, the limit is the GEM distribution with parameter 9/b. 
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1. Introduction 



We want to study and give some properties about several birth-death and immigration 
models where immigration is structured. All individuals behave independently of the others, 

f^ their lifetimes are i.i.d. but non-necessarily exponential and each individual gives birth at 

constant rate b during her life. We will consider the supercritical case i.e., the mean number 
of children of an individual is greater than 1. In the absence of immigration, if X{t) denotes 
the number of extant individuals at time t, the process (X(t),t > 0) is a particular case of 

^vq Crump- Mode- Jagers (or CMJ) processes [12] , also called general branching processes. Here, 

X is a binary (births arrive singly) and homogeneous (constant birth rate) CMJ process. 
Now, we assume that at each arrival time of a Poisson process with rate 9, a new individual 
enters the population and starts a new population independently of the previously arrived 

l/~) ones. This immigration model extends to general lifetimes the mainland-island model of S. 

i^ Karlin and J. McGregor p3j. In that case, the total population process X is a linear birth- 

and-death process with immigration. For more properties about this process, see [23] or [26] 
and references therein. In the context of ecology [1] , this model can be used as a null model 
of species diversity, in the framework of the neutral theory of biodiversity [11] . 

We first give the asymptotic behavior of the process (I(t),t > 0) representing the total 

S^ number of extant individuals on the island at time t. Specifically, there exists r/ > (the 

Malthusian parameter associated with the branching process X) such that e~ r,t I(t) converges 
almost surely. S. Tavare [22J proved this result in the case of a linear birth process with 
immigration. The case of general CMJ processes was treated by P. Jagers [12] under the 
hypothesis that the variance of the number of children per individual is finite. We manage 
to relax this assumption in the case of homogeneous (binary) CMJ-processes thanks to spine 
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decomposition of splitting trees [9] [16] [T7] which are the genealogical trees generated by 
those branching processes. In passing, we obtain technical results on the log-integrability of 
sup t e-^X(t). 

Then, we consider models where individuals bear clonally inherited types. They intend 
to model a metacommunity (or mainland) which delivers immigrants to the island as in 
the theory of island biogeography |19j . However, we made specific assumptions about the 
spectrum of abundances in the metacommunity. In Model I, there is a discrete spectrum with 
zero macroscopic relative abundances: when an immigrant enters the population, it is each 
time of a new type. In Model II, we consider a discrete spectrum with nonzero macroscopic 
relative abundances: the type of each new immigrant is chosen according to some probability 
(Pi, i > !)• In Model III, we consider a continuous spectrum of possible types but to enable a 
type to be chosen several times from the metacommunity, we change the immigration model: 
at each immigration time, an individual belonging to a species with abundance in (x, x + dx) 

is chosen with probability x g dx (where / is a positive function representing the abundance 
density and such that 9 := L xf(x)dx < oo) and it starts an immigration process with 

immigration rate x. The particular case of abundance density f(x) = ^—— appears in many 
papers. I. Volkov et al. [23] and G. Watterson [25] consider it as a continuous equivalent of 
the logarithmic series distribution proposed by R. Fisher et al. [8J as a species abundance 
distribution. In this particular case, species with small abundances are often drawn but they 
will have a small immigration rate. 

In the three models, we get results for the abundances P%, P2, ■ ■ ■ of different types as time t 
goes to infinity: the vector (Pi, P2, . . . ) rescaled by the total population size converges almost 
surely. More precisely, in Model I which is an extension of S. Tavare's result [22] to general 
lifetimes, we consider the abundances of the surviving families ranked by decreasing ages and 
the limit follows a GEM distribution with parameter 6/b. This distribution appears in other 
contexts: P. Donnelly and S. Tavare [5 J proved that for a sample of size n whose genealogy is 
described by a Kingman coalescent with mutation rate 6, the frequencies of the oldest, second 
oldest, etc. alleles converge in distribution as n — > 00 to the GEM distribution with parameter 
9; S. Ethier [7] showed that it is also the distribution of the frequencies of the alleles ranked 
by decreasing ages in the stationary infinitely-many-neutral-alleles diffusion model. 

In a sense, the surviving families that we consider in our immigration model are "large" 
families because their abundances are of the same order as the population size. A. Lambert 
[18j considered "small" families: he gave the joint law of the number of species containing k 
individuals, k = 1,2, . . . 

In Section[2j we describe the models we consider and state results we prove in other sections. 
Section [3] is devoted to proving a result about the process X, Section [4] to proving a property 
of the immigration process (I(t),t > 0) while in Section [5l we prove theorems concerning the 
relative abundances of types in the three models. 



2. Preliminaries and statement of results 
We first define splitting trees which are random trees satisfying: 

• individuals behave independently from one another and have i.i.d. lifetime durations, 

• conditional on her birthdate a and her lifespan £, each individual reproduces according 
to a Poisson point process on (a, a + £) with intensity b, 

• births arrive singly. 
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We denote the common distribution of lifespan £ by A(-)/6 where A is a positive measure on 
(0, oo) with mass b called the lifespan measure [T?] . 

The total population process (X(t),t > 0) belongs to a large class of processes called 
Crump-Mode-Jagers or CMJ processes. In these processes, also called general branching 
processes [121 ch.6], a typical individual reproduces at ages according to a random point 
process £ on [0, oo) (denote by /i := E[£] its intensity measure) and it is alive during a random 
time £. Then, the CMJ-process is defined as 

x ( f ) = ^2 *{«*<t<:a»+C*}> * - ° 

x 

where for any individual x, a x is her birth time and Q x is her lifespan. In this work, the 
process X is a homogeneous (constant birth rate) and binary CMJ-process and we get 



fx{dx) = dx I A(dr). 

J[x,oo) 

We assume that the mean number of children per individual m := J( 0oo \ rA(dr) is greater 
than 1 (supercritical case). 

For A > 0, define ip(X) := A— f, Jl — e~)A(dr). The function ip is convex, differentiable 

on (0, oo), ip(0 + ) = and ip'(0 + ) = 1 — L rA(dr) < 0. Then there exists a unique positive 
real number r\ such that ip{r]) = 0. It is seen by direct computation that this real number is 
a Malthusian parameter [12], p. 10], i.e. it is the finite positive solution of J °° e~ vr fj,(dr) = 1 



and is such that X(t) grows like e vt on the survival event (see forthcoming Proposition 2.1). 
From now on, we define 

c := tl)'(rf) 

which is positive because ip is convex. 

Another branching process appears in splitting trees: if we denote by Z n the number of 
individuals belonging to generation n of the tree, then {Z n ,n > 0) is a Bienayme-Galton- 
Watson process started at 1 with offspring generating function 



f(s) := [ 6- 1 A(dr)e- fer(1 - s) < s < 1. 

J(Q,oo) 



'(0,oo) 

To get results about splitting trees and CMJ-processes, A. Lambert [161 E] used tree 
contour techniques. He proved that the contour process Y of a splitting tree was a spectrally 
positive (i.e. with no negative jumps) Levy process whose Laplace exponent is tp. Lambert 
obtained result about the law of the population in a splitting tree alive at time t. If ¥ x denotes 
the law of the process (X(t),t > 0) conditioned to start with a single ancestor living x units 
of time, 

(1) F x (X(t) = 0) = W(t-x)/W(t) 

and conditional on being nonzero, X(t) has a geometric distribution with success probability 
l/W(t) i.e. fornGN*, 



w(t) J V w(t)J w(t) 
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where W is the scale function [2, ch.VII] associated with Y: this is the unique absolutely 
continuous increasing function W : [0, oo] — > [0, oo] satisfying 

f°° 1 

(3) / e~ Xx W{x)dx = —— X>7]. 

Jo WW 

The two-sided exit problem can be solved thanks to this scale function: 

(4) F(T <T (a ^ oo) \Y = x)=W(a-x)/W(a), < x < a 

where for a Borel set B of R, T B = inf{£ > 0, Y t £ B}. 

We now give some properties, including asymptotic behavior, about the CMJ-process X. 

Proposition 2.1. We denote by Ext the event \ lim X(t) = 0>. 

(5) P(Ext) = l-n/b 
and conditional on Ext c , 

(6) e-^Xit) -^4 E 



where E is an exponential random variable with parameter c. 
(ii) If, for x > 0, log + x := log x V 0, 



(7) E 



. 2 

log+sup(e- , ''X(t)) 
t>o 



Ext c 



< oo 



We remind that c = tp'{n) and then depends on the measure A(-). The proof of the last as- 
sertion requires involved arguments using spine decomposition of splitting trees. Proposition 
2.1 which will be proved in Sectional is known in a particular case: if the lifetime A(-)/6 



has an exponential density with parameter d, (X(t),t > 0) is a Markovian birth and death 
process with birth rate b and death rate d < b. In that case, r\ = b — d, c = 1 — d/b = P(Ext c ) 
and integrability of sup i>0 (e^ r ' i X(t)) stems from Doob's maximal inequality. 

We now define the immigration model: let 8 be a positive number and = Tq < T± < 
Ti < • • • be the points of a Poisson process of rate 9. At each time T«, we assume that a new 
individual immigrates and starts a new population whose size evolves like X, independently 
of the other populations. That is, if for i > 1 we call (Z l (t),t > 0) the ith oldest family (the 
family which was started at Tj) then Z l (t) = Xi(t — Tj)lr t >^.i where Xi, X2, . . . are copies of 
X and (Xi,i > 1) and (Tj, i > 1) are independent. This immigration model is a generalization 
of S. Karlin and J. McGregor's model [H] in the case of general lifetimes. 

For i > 1, define (Z^ l '(t),t > 0) as the ith oldest family among the surviving populations 
and T"W its birthdate. In particular, by thinning of Poisson point process, (T^ l ',i > 1) is a 
Poisson point process with parameter 6rj/b thanks to (pi). 

We are now interested in the joint behavior of the surviving families (Z^ l '(t),t > 0) for 
i > 1: 



e 
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Figure 1. splitting trees with immigration. The vertical axis is time, hori- 
zontal axis shows filiation. At time ti, three populations are extant. 



As in (pi), denote by Ei := lim^oo e~ vt Xu\(t) for i > 1. Thus E\, E2, ■ ■ ■ are i.i.d. exponential 
r.v. with parameter c. Moreover, the sequences (Ei,i > 1) and (T^\i > 1) are independent. 
It follows that e~ vt Z^ 1 ' (t) — > e _r? Ei a.s. as t — > 00. We record this in the following 

Proposition 2.2. 

e^\zW(t),Z^(t),...) -► (e-^^.e-^ ^,...) «... 

t— >oo 

where the Ei 's are independent copies of E and independent of the T^ 1 ' 's. 
For t > 0, let I(t) be the size of the total population at time t 

i(t) = ^z l (t). 

i>l 

The process (I(t),t > 0) is a non-Markovian continuous-time branching process with immi- 
gration. 

Theorem 2.3. (i) For t positive, I it) has a negative binomial distribution with parameters 
1 — H^(i) -1 and 9/b. i.e. for s £ [0, 1], its generating function is 

e/b 



G t (s):=E 



J(t) 



wity 1 



(ii) We have 



l-sil-Wit)- 1 ) 
I := lim e-^I{t) = V e^ W £i a.s. 

t — Inn £- — ^ 



»>1 



and 7 /tas a Gamma distribution T(6/b,c) i.e the density of I with respect to Lebesgue 
measure is 

„9/b x 9/b— lg— ex 



g{x) 



T(9/b) 



x > 0. 
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The result (i) is a generalization of a result by D.G. Kendall [15J which was the particular 
Markovian case of a birth, death and immigration process. The proof we give in Section [4] 
uses equations ([I]) and ^ about the law of X(t). 

There exist other proofs of the almost sure convergence in (ii), but they require stronger 
assumptions. For example, P. Jagers [12] gives a proof for the convergence of general branching 
processes with immigration under the hypothesis that the variance of the number of children 
per individual £(oo) is finite. In our case, this is only true if f, Q ^ r 2 A(dr) < oo. In the 

particular Markovian case described previously, the proof is also easier since (e~ vt X(t), t > 0) 
is a non-negative martingale [TJ p. Ill], (e~ r,t I(t),t > 0) is a non-negative submartingale and 
both converge a.s. In the proof we give in Section [4| the only assumption we use about the 
measure A is that its mass is finite. The proof is based on Proposition |2.1[ ii). 

In the following, we will consider different kinds of metacommunity where immigrants are 
chosen and will give results about abundances of surviving populations. In Model I, there is a 
discrete spectrum with zero macroscopic relative abundances: when a new family is initiated, 
it is of a type different from those of previous families. The following theorem yields the 
asymptotic behaviors of the fractions of the surviving subpopulations ranked by decreasing 
ages in the total population: 

Theorem 2.4 (Model I). 

limI(t)- 1 (Z^(t),Z^(t),...) = (P 1 ,P 2 ,...) a.s. 

t— >oo 

where the law of (P%, P 2 , ■ ■ ■ ) is a GEM distribution with parameter 9/b. In other words, for 
i> 1 

i-l 

Pi^BiJUl-Bj) 
i=i 

and (Bi)i>i is a sequence of i.i.d. random variables with law Beta(l,9/b) whose density with 
respect to Lebesgue measure is 

This result was proved by S. Tavare [22] in the case where A(dr) = 5oo(dr) (pure birth 
process); it is the exponential case defined previously with 6 = 1 and d = 0. His result is 
robust because we see that in our more general case, the limit distribution does not depend 
on the lifespan distribution but only on the immigration-to-birth ratio 9/b. In biogeography, 
a typical question is to recover data about population dynamics (immigration times, law of 
lifespan duration) from the observed diversity patterns. In this model, we see that there is 
a loss of information about the lifespan duration. However, the ratio 9/b can be estimated 



thanks to the species abundance distribution. We will prove Theorem |2.4| in Subsection 5.1 
In Model II, we consider a discrete spectrum with nonzero macroscopic relative abundances. 
Contrary to Model I where types were always new, they are now given a priori and types of 
immigrants are independently drawn according to some probability p = (pi,i > 1). When a 
population is initiated (i.e. at each time of the 6*-Poisson point process), it is of type i with 
probability pi > 0. 

Theorem 2.5 (Model II). For i > 1, denote by Ii{t) the number of individuals of type i at 
time t and set 9{ := -£*■. Then 

lim I(t)-\h(t)J 2 (t), ...) = (P[,P' 2 ,...) a.s. 

t— >oo 
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where for i > 1 

Pi = B'^l - B' 3 ) 

3=1 

and {B'i)i>i is a sequence of independent random variables such that 

B[ ~ Beta I 0j, - ^ pj 

\ j>i+l 



In particular, for i > 1, P[ has a Beta distribution B(9{, 9/b 



The proof of this theorem will be done in Subsection |5.2| In this model, the limit only 
depend on 9/b and the metacommunity spectrum (pi,i > 1). 

Remark 2.6. If the number of possible types n is finite, then 

n n t n O 

e^-Et- 1 - E«<=s 

and the joint density of (P{, . . . , P^) is 

T{9/b) 



72-1 \ 

n ^ _i ^>o} (i-^i — x n ^) e -- i \ xi „ 
i=i j 



1<1}- 



Hum) 

This is the joint density of a Dirichlet distribution Dir{9\, . . . , 9 n ). 

In Model III, we consider a continuous spectrum of possible types and we slightly modify the 
immigration process: when an individual arrives on the island, it starts a new population with 
an immigration rate proportional to its abundance on the metacommunity. More precisely, 
let II be a Poisson point process on M+ x IR + with intensity dt (8) xf(x)dx where / is a 
nonnegative function such that 9 := J °° xf(x)dx is finite. Then, write II := ((Tj,Aj),i > 1) 
where T\ < T2 < ■■■ are the times of a ^-linear Poisson point process and (Aj,i > 1) is 
a sequence of i.i.d random variables whose density is 9~ 1 xf(x)dx which is independent of 
(Ti,i > 1). At time Tj, a new population starts out and it evolves like the continuous-time 
branching process with immigration defined for the two previous models with an immigration 
rate Aj.The interpretation of this model is as follows: for x > 0, f(x) represents the density of 
species with abundance x in the metacommunity and at each immigration time, an individual 
of a species with abundance in (x, x + dx) is chosen with probability x q dx proportional to 
its abundance in the metacommunity. 

If (#*(£), t > 0) is the i-th oldest family, 

Z\t) = I i K {t-T i )l [t > Ti} , t>0 

where the i^.'s are independent copies of I&, which, conditional on A, evolves like the 
immigration process of the first two models with an immigration rate A. According to 
Theorem |2.3[ ii), we know that 

e-"*J A (t) — ► G a.s. 

t— s>oo 

where conditional on A, G ~ Gamma(A/6, c). We denote by F its distribution tail 



■:x 



T fM r°o p— ctj-x/b— l„x/b 



8 MATHIEU RICHARD 

Hence, we also have 
Proposition 2.7. For i > 1, 

e"^Z*(t) — » e-^Gi o.s. 

where (Gi,i > 1) is a sequence of i.i.d. r.v. with the same distribution as G and independent 
of{T h i>\). 

We again denote by I(t) the total population at time t: I(t) := \^ Z l (t) an d we obtain a 



result similar to Theorem 2.3 concerning the asymptotic behavior of I(t). 
Proposition 2.8. // J °° x 2 f(x)dx < oo we have 

e -*I(t) ±% y e -rfTi Gi 



i>\ 



and the Laplace transform of a := ^2i>i e Gj is 

E[e-]=exp(-^ °^(l-e-)^ 

Moreover, 

1 f°° 

E[a] = —— / x 2 f(x)dx < oo. 

We also have a result about abundances of different types 
Theorem 2.9 (Model III). We have 

'Z x {t) Z 2 (t) \ / o-i 0-2 



5 5 



o.s. 



J(i) ' J(i) ' " ' ) t^oo U',T 
where (o~i,i > 1) are £/ie points of a non-homogeneous Poisson point process on (0, oo) u>i£/i 
intensity measure ~-^p-dy and a = Yli>i a i- 



The proofs of these two results will be done in Subsection 5.3 Notice that in this model, the 
limit only depends on the lifespan measure via the Malthusian parameter n and the constant 
c. 

3. Proof of Proposition 12.11 

3.1. Some useful, technical lemmas. Thereafter, we state some lemmas that will be useful 
in subsequent proofs. 

Lemma 3.1. Let Y\,Y2, ... be a sequence of i.i.d random variables with finite expectation. 
Then, if S := sup n>1 (± Y!}=i Y i) , 



E 



(log + S)* 



< oo, k > 1. 



Proof. According to Kallenberg [T3l p. 184], for r > 0, 

rF(S> 2r) <E[Y 1 ;Y 1 > r\. 

Hence, choosing r = e s /2, we have for s > 0, 

P(log + S > s) < 2e~ s E[>i] 
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and 



poo poo 

E[(log+ S) k ] = / ks k - 1 F(log + S > s)ds < 2E[yi]fc / s k ~ 1 e~ s ds < 
Jo Jo 



oo. 



This completes the proof. 



□ 



Lemma 3.2. Let A be a homogeneous Poisson process with parameter p. If S := sup t>0 (At/t), 
then for a > 0, 



a 



In particular, 



\/k > 1, E 



(io g + sy 



< oo. 



Proof. If o < p, since lim^oo Aj/t = p, P(S' > a) = 1. Let now a be a real number greater 
than p. Then 

F(5 < a) = P(Vt > 0, at - A t > 0). 

According to Bertoin [2, chap.VII], since (at — A t ,t > 0) is a Levy process with no positive 
jumps and with Laplace exponent (f>(\) = Xa — p(\ — e _A ), we have as in (|4| 

where H is the scale function associated with (at — At, t > 0) . We compute H(0) and H(oo) 
using Tauberian theorems [2, p. 10]: 

• 4>(X) ~ X(a — p) then H (x) ~ (a — p)^ 1 



c/)(X) ~ Xa then H(x) ~ a 1 

oo o 



Hence, 



Then, 



^ \ a ~P 1 Z 3 

< a) = = 1 , a > p. 

a a 



poo 

E[(log+ S) k ] = / A;r fe - 1 P(log+ 5 > r)dr 
Jo 



< 



log+p 



kr k x dr + 



r k— ] P r 

fcr — dr < oo 



□ 



'0 -'log+p 

and the proof is completed. 

3.2. Proof of Proposition 2.1| (i). A. Lambert proved in [16] that P(Ext) = 1 — r]/b and in 
|17j that, conditional on Ext c , 

e- vt X(t) C 



E 



t— >oo 



where E is an exponential random variable with parameter c. To obtain a.s. convergence, 
we use |20(, Thm 5.4] where O. Nerman gives sufficient conditions for convergence of CM J 
processes to hold almost surely. Here, the two conditions of his theorem are satisfied. Indeed, 
the second one holds if there exists on [0, oo) an integrable, bounded, non-increasing positive 
function h such that 



E 



sup 

t>0 



h(t) 



< oo 
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where we recall that £ is the lifespan duration of a typical individual in the CMJ-process X. 
Then, choosing h(t) = e _r?t , this condition is trivially satisfied. 

The first one holds if there exists a non-increasing Lebesgue integrable positive function g 
such that 

(8) / -j-re^nidt) < oo. 

Jo 9\t) 

Taking g(t) = e~" 1 with 7] > j3 > and recalling that ji(dt) = f,. sA(dr)dt, we have 



f 

JO 



(*|t») 

— T e-"V(d*) ^ / e (/3 -" )4 / A(dr)di - / A(dr) /"" e^-^'dt 

5(*J JO Jt J(0,oo) Jo 



<c[ A(dr) = Cb<oo 

and condition ffify is fulfilled. 
3.3. Proof of Proposition |2Tl](ii) . We will say that a process G satisfies condition (C) if 



E 



log + sup e-^Gt] 
t>o J 



< oo 



and our aim is to prove that, conditional on non-extinction, the homogeneous CMJ-process 
(X(t),t > 0) satisfies it. 

According to Theorem 4.4.1.1 in [TS], conditional on non-extinction of (X(t),t > 0), 

x{t) = x? + xt + x? 

where 

• Xt° is the number of individuals alive at time t and whose descendance is infinite. In 
particular, (X 4 °°,t > 0) is a Yule process with rate rj. 

• Xf is the number of individuals alive at time t descending from trees grafted on the 
right hand side of the Yule tree (right refers to the order of the contour of the planar 
splitting tree) 

Nt 
Xt:=Y,Xi(t-fi) 

i=l 

where 

~ (Xi,i > 1) is a sequence of i.i.d. splitting trees conditional on extinction and 

independent of X°° . We know that such trees have the same distribution as 

subcritical splitting trees with lifespan measure A(dr) = e _r?r A(dr) (cf. |16j). 
— Conditionally on (X^° ,t > 0), (Nt, t > 0) is an non-homogeneous Poisson process 

with mean measure {b — rj)X^°dt and independent of (Xj)j. We denote its arrival 

times by Ti , T2 ; • • • 

• X® is the number of individuals alive at time t descending from trees grafted on the 
left hand side of the Yule tree (left also refers to the contour order) . 

More specifically, let (A, R) be a couple of random variables with joint law given 
by 

(9) P(A + R e dz, Redr) = e~^drA(dz), < r < z 
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and let ((Aij, Rij),i > 0,j > 1) be i.i.d random variables distributed as (A,R). We 
consider the arrival times 

T id = Ti + A iyl + Ai j2 + ■■■ + A itj , i > 0, j > 1 

where = tq < t\ < t 2 < ■ ■ ■ are the splitting times of the Yule tree, that is, on each 
new infinite branch, we start a new j4-renewal process independent of the others. We 
define for t > 0, 



E- 

i,3 



where (Xij,i > 0,j > 1) is a sequence of i.i.d. splitting trees independent of X°°, 
conditioned on extinction and such that the unique ancestor of JQj has lifetime Rij- 
We denote by N t := #{(i, j), Tij < t} the number of graft times before t. 



▲ I 



R-2 



Tl 



Ri 



^ 



Ti 



Rs 



a 



To, 2 T2 



To 



Ti, 



, r o,i 



Figure 2. Spine decomposition of a splitting tree. In bold, the Yule tree 
X°° on which we graft on the left (at times To.i, . . . ) the trees conditioned on 
extinction whose ancestors have lifetime durations distributed as R and on the 
right (at time Ti, . . .) the trees conditioned on extinction. 



In the following, we will prove that (Xt,t > 0) satisfies condition (C). To do this, using 
Minkowski inequality and the inequality 

Vx, y > 0, log + (x + y) < log + x + log + y + log 2, 

we only have to check that the three processes X°°, X 9 and X d satisfy (C). 

3.3.1. Proof of condition (C) for X°° . 

Since (X^°,t > 0) is a 77- Yule process, (e~' nt X£°,t > 0) is a non-negative martingale [U p. Ill] 

and so by Doob's inequality [21J, 



(10) 



E 



sup(e-" t Xr) 2 
.t>o 



<4supE[(e-' ?< X t 00 ) 2 ] 
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Moreover, E[(X°°) 2 ] = 2 (e 2 ^ - e^) (again in |lj) and 



E[(e" r ' t X 



oo\2i 



-2 V t (2 V t _ „nt 



(e 2r " - e*) 



2e- zr ' 1 (e zr i l - e T ' 1 ) — ► 2. 

t— >oo 



Hence, the supremum in the right hand side of (|10|) is finite and 

2" 



E 



snp(e-"%°°; 

4>0 



< oo. 



Prom now on, we will set M := sup i>0 e ^X^ . Since E[M 2 ] < oo, (C) is trivially satisfied 
by X°°. 

3.3.2. Proof of condition (C) for X . 
We recall that 

Nt 



xt = Y J x i {t-f i \ 



i=l 



Denote by Yi the total progeny of the conditioned splitting tree Xi, that is, the total number of 
descendants of the ancestor plus one. Then, a.s for all t > and i > 1, we have Xi(t— Tj) < Yi 
and 

Nt 
Xf <^Yi a.s. t > 0. 



i=l 



Hence, almost surely for all t, 



e^Xf < e~* 



Nt / N t 



i=l 



7V t 



t=l 



and, thanks to Minkowski's inequality, 



E 



log+supfe"^ 
t>o \ 



1/2 



< E 



log + sup ( e~ nt N t 



rn 



+E 



1/2 



N, 



log+ «>» (sS 15 



1/2 



We first consider the second term in the right hand side of (11): 

2n 



E 



log + sup -^ V Yi 
t>o \N t jri 



< E 



':°^M Y )) 



since (Nt,t > 0) is integer-valued. Thanks to Lemma 3.1, this term is finite because E[Yi] is 
finite. Indeed, Y\ is the total progeny of a subcritical branching process and it is known [10| 
that its mean is finite. 
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We are now interested in the first term in the r.h.s. of (11). We work conditionally 
on X°° =: (f(t),t > 0). Since we have e _r? ' L f(s)ds < M using the supremum M of 

rt 
e -*N t = e" 



-</' / f {s ) ds ^— < M^^- 



fof(s)ds fof(s)di 



Moreover, 



N t} t>0 



(j) 



Jo /(*)<**' 



°) 



where N' is a homogeneous Poisson process with parameter b — rj. 
Hence, using Minkowski's inequality, 



E 



log + sup (e'^Nt 
t>o \ 



1/2 



< log + M + E 



log + M + E 



i -i- I fof( s ) ds 

log^sup I — y 



t>0 \ Jo /( s ) di 



1/2 



'l + f^t 

log^ sup — 



1/2 



and the second term of the r.h.s. is finite using Lemma 3.2 
Hence, (iVj, t > 0) satisfies (C) since E[M 2 ] < oo and X rf as well, which ends this paragraph. 

3.3.3. Proof of condition (C) for X 9 . 
We have 



x t =j2Mt-Ti) 



i=l 



As in the previous section, 



/ 1 M 

\ * i=l 



a.s. 



where % is the total progeny of the conditioned CMJ-process (Xi(i),t > 0). 
Hence, 



E 



log+supfe-^Xf)^) 



1/2 



< E 



(12) 



log+supfe-^iVtH 

: og+ s(S"')) 



1/2 



We first prove that the second term in the r.h.s. is finite using Lemma |3.1| We only have 
to check that E[Yi] is finite. We recall that Y\ is the total progeny of a splitting tree whose 
ancestor has random lifespan R\ and conditioned on extinction. Conditioning on R±, it is also 
the total progeny of a subcritical Bienayme-Galton- Watson process starting from a Poisson 
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random variable with mean Ri. Hence, E[Yi] is finite if E[i?i] is finite. As a consequence of 
Q, we have that F(Ri G dr) = e~^ r J r °° K{dz)dr and 

E[i2i] = /" A(efe) [ re-^dr 
J(a,oo) Jo 

A(dz)- +/ A(dz) 2 

(0,oo) V J(0,oo) V 

77 T) T) 



We are now interested in the first term of the r.h.s. of (12). We need to make calculations 



on iVj which is the total number of times of graftings Tij less than or equal to t. Recall that 
for i > 0, j > 1, Tij = t% + Ajj where Aij := A^i + • • • + Ajj and that Tj is the birth time 
of individual i and that Nt is the sum of the numbers of graftings before t on each of the 
Xf° branches. For i > 0, denote by a|, a^, ■ • • , the birth times of the daughters of individual 



i and a Q 



t%. For k > 1, denote by f£ 



°4 ~~ a fe-i the interbirth times. In particular, 
(fl,i > 0, k > 1) are i.i.d. exponential r.v. with parameter 77 since we consider a 77- Yule tree. 



Ai A 2 Oil A 3 A 4 oi 2 A 5 A 6 A 7 a 3 



Ai 



4i 



A 2 



A, 



As 



Figure 3. Construction of the renewal process (Aj,j > 1) by concatenation 
of the renewal processes Ak- 



We enlarge the probability space by redefining the renewal processes (Aij,j > 1) from a 
doubly indexed sequence of i.i.d. A-renewal processes (Ai t k,i > 0,/c > 1). We define the 
process (Aij,j > 1) recursively by concatenation of the Ai t k s as m Figure 3 To simplify 
notation, we only define (Aqj, j > 1) which will be denoted by (Aj,j > 1). Then, 



and 



A 1 := inf{t > ai\Ai n [t, +00) / 0}, C\ := #Ai n [0, ai] + 1, 



A,- := inf{t > 0|#Ai n [0,i] = j}, j = 1, . . . , Cy + 1. 



Moreover, for / > 1, if one knows C\ and A;, let r\ be the unique integer such that A ; belongs 
to ]a n ,a n+ i] and define 

A' l+1 := A[ + inf{< > a n+ i - A,|A n +i n [t, +00) / 0}, 
C /+ i := #A,+i n [0, a n+ i - A[] + 1, 
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in 



and 

A Cl +...+c l+ j ■= X + inf{t > 0|#A i+ i n [0,t] = j}, j = l,..., Q+i. 

Then, (Aj,j > 1) is a A-renewal process because we concatenated the independent renewal 
processes (Ak) stopped at the first renewal time after a given time. Indeed, one can see a 
renewal process as the range of a compound Poisson process whose jumps are distributed as 
A\. The renewal process stopped at the first point after t is then the range of a compound 
Poisson process stopped at the first hitting time T of [t, oo), which is a stopping time. Then 
(Aj,j > 1), as the range of the concatenation of independent compound Poisson processes 
killed at stopping times, is a compound Poisson process, by the strong Markov property. In 
conclusion, (Aj,j > 1) is a ^4-renewal process. 
According to previous computations, we have 



N t < 



E 

i>0i>l 



C?lK u <t} 



where C\ and r); are the analogous notations as C\ and r% for individual i. Then if 



Dl:=#A iik n[0,fl] + l k>l,i>0, 



since a 



n,i+i 



A ti <a 



l.i 



(13) 



n,i+i 

N t < 



a 



n. 



ft , then Cf < Dt and 



i>0,l>l i>0.k>l 



^k\a\<t} a - s ' 



Moreover, the random variables (D l k ,i > 0, k > 1) are independent and identically dis- 
tributed as 1 plus the value of a ^4-renewal process at an independent exponential time E 
with parameter 77 

D := sup{j > 1, A x + • • • + Aj < E} + 1. 

1 terms. Indeed, each individual of the Yule tree 



The sum of the r.h.s. of (13) has 2X t °° 



contributes to 1 + rtj terms in the sum where n% is the number of daughters of individual i 

rii terms in the sum and Y^i n % ^ the number of 

1. 



born before t. Then, there is X^ 



E 

j=0 



descendants of the individual born before t which equals X\ 
Hence, using Minkowski's inequality, 



E 



log^ 



sup 

4>0 



-tfjVi 



1/2 



< E 



log H 



sup [e 



-T)t( 



{2X? - 1)) 



1/2 



(14) 



+E 



log sup 

n>0 




1/2 



where (D{,i > 1) are i.i.d. r.v. distributed as D. The first term of (14) is smaller than 

if E[Di] < 00. 



3.1 



E[(log + (2M)) 2 ] 1 ' 2 and the second term is finite by another use of Lemma 
Denote by (Ct,t > 0) the renewal process whose arrival times are distributed as A. A little 
calculation from ft9h gives us that EL4] = ^—— > 0. Thus, according to Theorem 2.3 of 
Chapter 5 in [B], we have 

lim E[01 = ^ 

i->oo t 771—1 
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and so there exists k > such that E[Ct] < nt for t > 0. Then, 

/•oo />oo 

E[£>i] = 1 + E[Cy = 1 + / ^e^diE^t] <1 + k rje^dt = 1 + - < oo. 

7o 7o »7 



Finally, the r.h.s. of (14) and (12) are finite and X 9 satisfies condition (C). 

4. Proof of Theorem 12.31 

4.1. Some preliminary lemmas. We start with some properties about W the scale function 
associated with ifi and defined by pi). 

Lemma 4.1. (i) W(0) = 1, 

(ii) e'^Wit) — > c- 1 as t ^ oo, 
(iii) W is differentiable and W -k A = bW — W where * is convolution product. 



Proof. (i) We have 



OO J 1 

e- xt W(t)dt=— -r ~ - 

tp(A) A->oo A 



because ip(X) = A — b + f, Q , e Xr A(dr). Then, by a Tauberian theorem [2j plO], 
lim t ^ W(t) = l. 
(ii) (From [T7]) For A > 0, using a Taylor expansion and ip(rj) = 0, 

t/>(A + rj) ~ Xip'(rj) = Ac. 

A— >0 
Then 



/•OO 1 

/ W^e-^e-^dt ~ - 
7o A ^° A< 



A^O Ac 

and another Tauberian theorem entails that W(t)e~ rit converges to 1/c as t — > oo. 
(iii) We first compute the Laplace transform of W * A. Let A > r/ 

/•OO /•OO /• 

/ e~ xt W * A(t)dt = e~ xt W(t)dt e~ Xr A(dr) 

JO JO J(0,oo) 



Integrating by parts and using (i) and (ii), 



OO 



'IX 



e - xt W'(t)dt = e~ xt W(t) + A / e~ xt W{t)dt = -1 + 
o 



oc 



A 



o Jo " "' ' V'(A) 

and so the Laplace transform of bW — W is b/ip(X) + 1 — X/tp(X) which equals that of 
W -k A. This completes the proof. 

□ 

The following lemma deals with the convergence of random series: 

Lemma 4.2. Let (Q,i > 1) be a sequence of i.i.d. positive random variables such that 
E[log + Ci] is finite and let (tj,z > 1) be the arrival times of a Poisson point process with 
parameter p independent of (Q, i > 1). Then for any r > 0, the series YLi>i e ~ rn d converges 
a.s. 
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Proof. We have 
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(15) 



We use the following consequence of Borel-Cantelli's lemma: if £1,^2, ■■■ are i.i.d. non- 
negative random variables, 



1. sn 

lim sup — 

n— >oo Tl 



or 00 a.s. 



according to whether E[£i] is finite or not. We use it with £j = log + Q. Hence, since E[log + Ci] 
is finite, lim.j_-.oo log + Q/i = a.s. Moreover, by the strong law of large numbers, Ti/i converges 
almost surely to 1/p as i goes to infinity. Then, 



n iog + d 



> a.s. 



i— >oo p 



So the series in (15) converges a.s 



□ 



4.2. Proof of Theorem 2.3| (i). In order to find the law of I(t) the total population at 
time t, we use the fact that it is the sum of a Poissonian number of population sizes. More 
specifically, if we denote by Nt the number of populations at time t, (Nt,t > 0) is a Poisson 
process with parameter 9 and conditionally on {Nt = k}, the /c-tuple (Ti,...,Tfc) has the 
same distribution as (C^(i), • • • , ^(fc)) which is the reordered /c-tuple of k independent uniform 
random variables on [0,t]. Hence, conditionally on {Nt = k}, 

i( t ) ( i $>(*- u w) - IX* - °i) = XX^) 



4 = 1 



i=l 



i=l 



since all Ur^'s appear in the sum and the C/j's are independent from the Xj's. Hence, 
conditionally on {Nt = k}, I(t) has the same distribution as a sum of k i.i.d. r.v. with law 
X(U). Then, 



G t (a) =£>[.• 



J(t) 



N t = k 



W(N t = k) 



k>0 



(16) 



^[s Xl{Ul) 

k>0 

e~ et exp (tffE 



k (t0) k -et 
k\ ' 



We now compute the law of X(t) for t > and then we will compute the law of X\{U\) 
Using (fTl), Q and Lemma |4.1 [in), we have 



F(X(t) = 0) 



A(dr) 



F r (X(t) = 0) 
(0,00) 

1 



W(t - x) A(dr) 
.0,oo) ^(*) 6~ 



6W(t) 



W*A(<) = 1 



lU'(t) 
6W(t) 
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and for n G N*, 

F(X(t) = n) 



(0,00) 

1 

bW(t) 



¥ r (X(t) = n 

1 
1 



1 



1 



W(t) 

n— 1 



A(dr) 
b 

n-l 



6- W(t) _1 W*A(t) 



W(t)y bW(t) 2 ' 



We now compute P(Xi(C/i) = n) 

¥(X l (U 1 ) = 0) = 



If 1 If 1 W'(u) 

F(X 1 (u) = 0)du = l-~ ^fdu 



1 



logW(t) 
tb 



because W(0) = 1. For n > 0, 

1 * 



P(X 1 (C/i) = n) 



1 
bt 



F(Xi(u) = n)du 



1 



1 



1 



(1-u) 



n-l 



-du 



w(t)- 



n 



;i - i/^(*)) w 

btn 



n-l 



6VF(w) 2 



rfu 



We are now able to compute the generating function of X\(Ui). For s,t > 0, 



E 



,*i(tfi) 



1 



£ 

n>l 



S 

n 



1 



W(t) 



+ 1 



log l^(t) 



^ [log ( 1 -* ( .L-.i/Tn/)))+ logino 



1 



1 

6i 



bg(W(t) + a(l-W(t)) 



Finally for £, s > 0, according to (16) 
G t (s) = e- et exp(t9E 



*Xi(Ui) 



(W(t) + s(l-W(t)))- e/b . 



which is the p.g.f. of a negative binomial distribution with parameters 1 — W(t) and 9/b. 



4.3. Proof of Theorem |2.3| (ii). We first prove the almost sure convergence. Splitting I(t) 
between the surviving and the non-surviving populations, we have 

(17) e~^I(t) = J>"*zW(t) +'£ f e-*X i (t - T^ronExt, 

where for i > 1, Extj denotes the extinction of the process Xj. We will show that for each 
of these two terms, we can exchange summation and limit, so that in particular, the second 
term vanishes as t — ¥ 00. 

We first treat the second term of the r.h.s. of (|17|). We have 



Ct ■— 2_^Xi(t - Tj)l{ t > T .} nExti < 2^ l{ t > T .}yjl E xti a.s. 



i>\ 



i>l 
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where Yi is the total progeny of the i-th population which does not survive. Moreover, 
E [YilExtJ ^ JE[Yi]> Y\ is the total progeny of a subcritical Bienayme-Galton- Watson process 
so its mean is finite. Hence, since the r.h.s. in the previous equation is a compound Poisson 
process with finite mean, it grows linearly and e~ vt Ct vanishes as t — > oo. 

To exchange summation and limit in the first term of the r.h.s. of (17), we will use the 
dominated convergence theorem. By Proposition 2.2 we already know that e~ rit Z^ l \t) a.s. 
converges as t goes to infinity to e~ v Ei. Hence, it is sufficient to prove that 



(18) 
Since 

we have 



^sup(e- r ' 4 Z»(t)) <ooa.s 



j>i 



t>o 



supie-"*zW(t) 



e -^ W sup(e-^X (i) (t)) 

4>0 



j>i 



sup ( e 



i*zW(t) 



bTW 



»>1 



J, 



4.2 



where Jj := sup t>0 (e r,t X(^{t)) for i > 1 and J\, J2, ■ ■ ■ are i.i.d. Thus, using Lem mas 
this series a.s. converges if E[log + J\] is finite, which is checked thanks to Proposition |2. l[ ii) . 
Then we get (18) and using the dominated convergence theorem, 

I := lim e-^Ut) = V e^ E { a.s. 



i>\ 



In order to find the law of /, we compute its Laplace transform. For a > 0, using part(i) 
of this theorem, 



E 



-ae-^I(t) 



G t {e- 



-vt. 



-vt 



+ 1 



—ntt-nt 



) W(tj) 



-e/b 



and 



-v* 



+ 1 



-17* 



W(t) ~ 1 + ae~ vt W{t) — >■ 1 + - 

t— >oo t— >oo c 



using Lemma 4.1 (ii). Then, 



E[. 



-a/1 



a + c 



e/b 



which is the Laplace transform of a Gamma(0/6, c) random variable. 

5. Other proofs 



5.1. Proof for Model I. To prove Theorem 2.4, we will follow Tavare's proof |22j. We begin 
with a technical lemma which will be useful in the proof of this theorem and in forthcoming 
proofs. 

Lemma 5.1. Let (T%,i > 1) be the arrival times of a Poisson process with parameter p 
and d,^2,--- be i.i.d. r.v. independent of the Ti's. Denote by g the density of £1 with 
respect to Lebesgue measure and by F(v) := P(Ci > v) its distribution tail. Then for r > 0, 
(e~ rTi Ci,i > 1) are the points of an non-homogeneous Poisson point process on (0, 00) with 
intensity measure -—^-dv. 
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Proof. We first study the collection II = {(Tj,^),i > 1}. This is a Poisson point process on 
(0,oo) x (0,oo) with intensity measure pg{y)dtdy. Then, (e~ rTi Ci,i > 1) is a Poisson point 
process whose intensity measure is the image of pg(y)dtdy by (t, y) t-t e~ rt y. We now compute 
it. Let h be a non- negative mapping. Changing variables, we get 



h(e rt y)pg(y)dtdy 



r 


/■DC 

Jo 


h(v)dv 


'0 





du 
u 2 


P 

r 


/"OC 

Jo 


h{v) F{v) dv 

V 







and the proof is completed. 
We are now able to prove Theorem |2.4| By Proposition |2.2| and Theorem 2.3 

e-" ( (ZW(t),Z( 2 )(t),...) 



□ 



t->oo \ a <t 



a.s. 



i{t)-\zU{t),zW{t),..:, ( _ ////(/) 

where <7j := exp (— 77TW) E^ and <r := Ej>i a i- 

Moreover, the (<Tj)i>i are the points of a non-homogeneous Poisson point process on (0, 00) 

~ 9 e" cy 

with intensity measure dy thanks to Lemma 

b y 

because (T^')i>i is a Poisson process of rate 9r]/b and (Ei)i>i is an independent sequence of 

i.i.d. exponential variables with parameter c. 



5.1 



with p = 9i]/b, r = 7} and F(v) 



According to [3, p. 89], the Poisson point process (<7i)j>i satisfies a 
(actually a has a Gamma distribution) and the vector 

V a ' a ' "/ 

follows the GEM distribution with parameter 9/b and is independent of a. 



E 



i>\ 



Oi < OO 



5.2. Proof for Model II. We will prove Theorem [23} We recall that in Model II, immigrants 
are of type i with probability pi. Denote by N l (t) the number of immigrants of type i which 
arrived before time t. Then, (N l (t),t > 0) is a Poisson process with parameter 9pi and 
the processes (N l ,i > 1) are independent. Hence, Ii(t), hit), . . . are independent and their 
asymptotic behaviors are the same as I(t) in Theorem |2.3[ replacing 9 with 9pi. Then we 
have 



-*Ht) 



I; 



t— >oo 



a.s. 



i > 1 



where the Ij's are independent and I{ has a Gamma distribution T(9i,c) (recall that 9i 
9 Pl /b). 
Moreover, 

e-^I{t) — »■ / a.s. 



t— >oo 

where I ~ T(9/b, c). Therefore, for r > 1, 

lim J(t)- 1 (/ 1 (t),...,J r (t)) = 

i— ^00 

In order to investigate the law of this r-tuple, we prove that 
(19) I = Y^Ii a.s. 



a.s. 
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First, by Fatou's lemma, 



«>i 



and so 



Second, 



liminfe - ''* Vij(i) > Vliminf e^Idt) a.s. 



a.s. 



i>\ 



E 






Et 



i>l 



be 



E[I]. 



The last two equations yield (19). For 1 < i < r, we can write 

h 

I 



h 



h + • • • + l r + 1* 

where I* is independent of (Ij, 1 < % < r) and has a Gamma distribution T{6/b — 9 r ,c) with 
9 r := J2i=i @i- Hence, one can compute the joint density of the r-tuple (h/Ii ■ ■ ■ , Ir/I) as 
follows 

T(9/b) 



f{xi,...,X r ) 



wb-Ormum 



6*1—1 r —l(-\ _ 



Ju J 



\e/b-e r -i 



for xi, . . . , x r > satisfying x\ + - ■ - + x T < 1. This joint density is exactly that of (P{, . . . , P' r ) 
defined in the statement of the theorem. 



5.3. Proofs for Model III. We first prove the almost sure convergence in Proposition 2.8 



In order to do that, we use the same arguments as in the proof of Theorem |2.3[ ii): we will 
use the dominated convergence theorem for the sum 



e-" t /(t) = ^e-^/X s (t-T i )l {t > Ts} . 
As in a previous proof, this sum is bounded by 2, e ~ v " SU P ( e_r?< ^A;(*)) which, according to 



i>i 



Lemma 4.2, is a.s. finite if 



E 



t>o 



< oo. 



log+sup(e-^4(i)) 
*>o 

However, I^(t) = > X l (t — Tj)lr 4> j,| where conditionally on A, (Ti,i > 1) is a Poisson 

i>l 

process with parameter A. Hence, 

sup (e-^I^t)) <J2 e ~ vft SU P {e- vt X l (t)) =Y,e- vf *Ji 



*>o 



i>l 



t>0 



i>l 



where J±, J2, . . . is an i.i.d. sequence of random variables independent from Ti,T2, . . . dis- 
tributed as sup i>0 (e _,? *X(t)) where (X(t),t > 0) is a homogeneous CMJ-process. According 
to Proposition 2.1 'ii), we know that E[(log + J1) 2 ] < 00. 

We define for i > 1, q := e~ vTi Ji and ? := J2i>i ft an d we have to prove that E[log + c] 
is finite. To do that, we first work conditionally on A. According to Lemma 5.1 we know 
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that ($i,i > 1) are the points of a non-homogeneous Poisson process on (0, oo) with intensity 
measure —dv where L(v) := P(J > v). Then, using the inequality 

log + (x + y) < log + x + log + y + log 2, x, y > 0, 



we have 

(20) E[log + ft] <log2 + E 



log+^Ql^ 



<i} 



i>l 



+ E 



We first consider the second term of the r.h.s. 



(21) 



E 



log + J^Ql{ ft <i} 



«>i 



< E 



£«* 



%<i} 



«>i 



log+^ftl^x} 



1 AL(«) , A 
v —dv< —. 
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Then, we compute the third term of the r.h.s. of (20): if A := supj ft, 



E 



log+^ftl {ft 



>i} 



i>l 



<E[log+(^-#{i>l| Q >l})] 

< E [log+ A) + E [log + #{* > l|ft > 1}] 



Furthermore, the number of ft greater than 1 has a Poisson distribution with parameter 
Ji°° f^*;. Since E [ lo S + J ] < °°> 



/"OO /*00 /»C 

/ P(log+ J > s)ds = P( J > e s )ds = / 

Jo Jo Jl 



L(v) 



dv < oo. 



Then, 



A f°° L(v) 
(22) E [log+ #{i > l|ft > 1}] < E [#{i > l|ft >!}]=-/ -U-dv < CA 



V Ji v 



where C is a finite constant which does not depend on A. 
We now want to study A: 



/ roc a L( v ) \ 
F(A <x) = P(#{i > l|ft > x} = 0) = exp I - / y -!-dv J 



dv ) , x > 0. 



So that 



Then, 



7] X 



AL(v) 



rj v 



dv ) dx. 



E[log+ A] 



AL(x) 

logx exp 

77 x 



AL(v) 



x V v 
00 



dv dx 



(23) 



where C is a finite constant since E 



< - f°° logx^-dx = - I uL(e u yiu 

V Ji x r, J 

A f°° 

< — / uP(log+ J > u)du < C'A 

V Jo 

(log + J) < 00 according to Proposition 



with (21), (22) and (23) we have 



2.1 



'ii). Hence, 
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E[log + ?|A] <log2 + C"A. 
Then, E[log + ?] is finite because E[A] = 9~ l J °° x 2 f(x)dx < oo. and, using the dominated 
convergence theorem, e~ vt I(t) a.s. converges toward a = X]i>i e~ vTi Gi as t — > oo. 

We now compute the law of the limit a. We define cij := exp (— 7/Tj) Gi for i > 1. Then, 
using Lemma 5.1 , (<7j)i>i are the points of a non-homogeneous Poisson point process on (0, oo) 



with intensity measure -^-^-dy where F(y) = P(G7 > y). To compute the Laplace transform 



9 F(y) 
V V 

of a, we use the exponential formula for Poisson processes: for s > 0, 

E[e-]=e Xp (-^ °^M(l-e-)^). 
and to get the expectation of a, we differentiate the last displayed equation at 0: 

E[a] 
and 

E[G\: 

Hence, 



F(v)dv = -E[G1 



xf(x) — -dx < oo 
o be 



EH 



1 
rjbc 



x f(x)dx < oo 



which ends the proof of Proposition 2.8 



It remains to prove Theorem 2.9 that is to show that the vector (Z 1 (t), Z 2 (t), . . . )/I(t) a.s. 



■F(y). 



converges to a Poisson point process with intensity measure — —dy. It is straightforward 



using previous calculations, Propositions 2.7 and 2.8. 



Remark 5.2. Thanks to similar calculations as in Theorem, 2.3(i), we can compute the 
generating function of I(t) 



E 



Mt) 



exp 



du\ 9 



xf(x) 



dx 



, se[0,l]. 



io \ Jo (W(u)(l-s) + s) x/b )) 
and we can deduce the law of a in another way. 
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